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RESEARCH MEMORANDUM 

THEORETICAL AND EXPERIMENTAL INVESTIGATION OF THE 

EFFECT OF YAW ON HEAT TRANSFER TO CIRCULAR 

CYLINDERS IN HYPERSONIC FLOW 

By A. J. Eggers^ Jr.^ C. Frederick Hansen^ 
and Bernard E. Cunningham 


SUMMARY 


An approximate theory is developed for predicting the rate of heat 
transfer to the stagnation region of hlunt "bodies in hypersonic flight. 
Attention is focused on the case where wall temperature is small compared 
to stagnation temperature . The effect of yaw on heat transfer to a cylin- 
drical stagnation region is treated at some lengthy and it is predicted 
that large yaw should cause sizable reductions in heat-transfer rate. 

Experiments were conducted in a so-called hypersonic gun tunnel to 
obtain a preliminary check on this theoretical prediction. These experi- 
ments consisted of measuring the rate of heat transfer to circular cylin- 
ders in a hypersonic air stream of nominal Mach number^ 9*8^ and nominal 
stagnation temperature, 2200*^ R. Experiment tended to confirm theory, 
showing, for example, that a 60 -percent reduction in average heat-transfer 
rate can be obtained by yawing a circular cylinder 70 ^ to the air stream. 


INTRODUCTION 


It has been suggested (see refs. 1 and 2) that blunting or rounding 
the leading edges of wings and bodies might substantially alleviate aero- 
dynamic heating of these regions in hypersonic flight. There is, of 
course, the added advantage that round leading edges are structurally 
more practical than sharp leading edges, especially when the problem of 
absorbing heat is considered. Another consequence of round leading edges 
may be increased pressure drag. In the case of ballistic vehicles, this 
consequence is often an advantage (see ref. l) . In the case of glide 
vehicles, however, or more generally any vehicles required to operate for 
sustained periods in more or less level hypersonic flight, Increased drag 
would most likely be viewed as a disadvantage. 
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Now, to "be sure, rounding or "blunting the nose of a hody does not 
always increase drag. Indeed, small amounts of blunting may reduce the 
drag of a hody (see, e.g., refs. 3 and i)-) . The same, however, cannot he 
said for blunting the leading edge of a wing. Even small blunting causes 
a sizable increase in drag. It is natural, then, to look for methods of 
minimizing this drag penalty, and the possibility of yawing or sweeping 
the leading edge comes to mind. Impact pressures should be, according to 
simple-sweep theory, decreased in proportion to the cosine squared of the 
angle of sweep; hence, as is intuitively obvious, large sweep sho\ald sub- 
stantially reduce the drag penalty due to blunting. In view of this pos- 
sibility it is important to inquire of the effect of yaw or sweep on heat 
transfer to a blunt leading edge. 

A theoretical study of this effect was therefore undertaken, and the 
predictions of the theory were checked against the results of experiments 
carried out in a hot hypersonic air stream. The purpose of this report 
is to describe the problems and resiUts of both the theoretical and experi- 
mental research. 


SYMBOLS 


A,B,C, 

D,E,F, 

G, • . . 

. integration 

constants 

Cb 

specific heat of body material, ft-lb/slug °R 

Cp 

specific heat at constant pressure, ft-lb/slug °R 

I>b 

2 R-| 3 , twice the radius of curvature of the body at the stagnation 
point , ft 

fl(t) 

pdt 

parameter of time, sec 

p(t) - p(o) 

f2(t) 


parameter of time, sec 

p(t) - p(o) 


H reservoir pressure, Ib/ft (unless otherwise specified) 

h specific enthalpy, ft-lb/slug 

k coefficient of thermal conductivity, ft-lb/ft-sec °R 

K ratio of static pressure to pitot pressure, dimensionless 
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M 

n 

Nu 

P 

Pt 

Pr 

q 

R 

Rb 

Rs 

Re 

r,0, cp 
T 
Tb 
To 

To(^) 

Tr 

Tt 

Uoo 

u,v,w 

u,v 


Mach mimher^ dimensionless 

exponent of temperature in thermal conductivity and viscosity 
functions (see eqs. (37) stnd ( 38 )), dimensionless 

Nusselt number based on a length D]^ and stagnation temperature 
conditions^ dimensionless 

static pressure, Ib/ft^ (unless otherwise specified) 
pitot pressure, Ib/ft^ (unless otherwise specified) 

Prandtl number, dimensionless 
dynamic pressure, Ib/ft^ 
gas constant, ft-lb/slug °R 

radius of curvature of body at the stagnation point, ft 

radius of curvature of the shock wave at the stagnation stream- 
line, ft 

Reynolds number, based on twice the radius of curvature of the 
body at the stagnation point, dimensionless 

spherical coordinates, feet, degrees, and degrees, respectively 

static temperature, °R 

temperature of the body , °R 

temperature at the interface, x = 0, with body at zero yaw, °R 

temperature at the interface, x = 0, with body at angle of yaw A 
°R 

recovery temperature , °R 
stagnation temperature, °R 
stream velocity, ft/sec 

velocity components in the x,y, and z directions, respectively 
ft/sec 

velocity components in the x and r directions, respectively, 
ft/sec 
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x,r 

6 

e 

7 


r 

7 

P 

P- 


a 

M- 

M-o 

^o(^) 

T 


S 

■b 

o 


00 


Cartesian coordinates, ft 
cylindrical coordinates, ft 
flow deflection angle, deg 

dimensionless coordinate, — or — ^ 

ratio of specific heat at constant pressure to specific heat at 
constant volume, dimensionless 

7 m 

7-1 

, a function of 7 , dimensionless 

angle of yav, deg 
density, slugs/cu ft 

density of the hody material, slugs/cu ft 
T 

k dT, a function of the coefficient of thermal conductivity 
and of temperature, ft-lh/ft-sec (unless otherwise specified) 
acute angle of shock wave relative to stream velocity vector, deg 
coefficient of viscosity, slugs/ft sec 

coefficient of viscosity at temperature Tq, slugs/ft sec 
coefficient of viscosity at temperature Tq(?\), slugs/ft sec 
time constant, sec 




Subscripts 

conditions just behind shock wave on the stagnation streamline 

conditions at the stagnation point of the body 

conditions at the interface between regions 1 and 2 on the stag- 
nation streamline (see sketch (a)) 

conditions in the free stream 
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Superscripts 


first derivative with respect to the x coordinate 
second derivative with respect to the x coordinate 


THEORY 

General Equations in Cartesian Coordinates 


The analysis proceeds from the equations of momentum, continuity, 
energy, and state for continuum fluid flow. The x, y, and z momentum 
equations are, respectively. 
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^5? 





_ 2 i_ 

r„^ 



3 ^z 


\dx dy dz / _ 


Sz 


Sw 




) dx 







(3) 


The continuity equation is 


dt 


+ ^(pu) + ^(pv) + :^(pw) = 0 
dx dy dz 


(M 


and the energy equation is 





^h 

bt 


)- 


( 




+ u 


dx 


+ V 





2 




+ 



du dv \ 
vSy dxy 


+ 



dv dwA 
.Sz dyy' 


2 ( ^ ^ ^ ) 
3 \dx dy dz J 


(5) 


while the equation of state is taken in the form 


p = p(p,T) 


( 6 ) 
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Derivations of the momentum and energy eq_uations are given in numerous 
sources (see^ e.g.^ refs. 5^ 6^ and 7). Note that the coefficients of 
viscosity and thermal conductivity^ and the heat capacity have been 
treated as variables. It is intended that by so doing a more accurate 
solution will be obtained for hypersonic flows with their characteristi- 
cally large temperature and pressure gradients. 

Let us now consider the particular flows of interest in this paper ^ 
namely ; those in the region of a stagnation point. 


Model of Flow and Method of Analysis 


It is instructive in setting up the model to consider the qualitative 
aspects of temperature and velocity variations in the flow along the stag- 
nation streeanline . Restricting the analysis to steady hypersonic flow, 
that is L^sin 6 >>1, we will assume that the surface temperatirre is low 
compared to the stagnation temperature of the air. This assumption seems 
quite reasonable since practical surface materials will probably be 
destroyed if surface temperat-ures are allowed to approach stagnation 
temperature. It will be assumed further that the Reynolds number of the 
flow is large enough so that heat conduction and viscous shearing in the 
shock process is distinct and separate from the corresponding phenomena 
occurring in the boundary layer adjacent to the surface of the body. 
Accordingly, temperature and velocity should vary along the stagnation 
streamline similar to the manner shown in sketch (a). 
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There is an abrupt and large increase in temperature and decrease in 
velocity of the air as it passes through the bow shock. Proceeding from 
the shock in the direction of the body, temperature continues to increase 
slowly while the velocity decreases slowly towards zero. Near the surface 
of the body, the air tempera t\ire ceases to increase and, in fact, begins 
to fall off steeply in the direction of the body temperature. The veloc- 
ity of the flow must, of co\irse, be close to zero in this region. 

On the basis of these observations the following simplified model is 
proposed and employed throughout this study of heat transfer in a stagna- 
tion region. 


Since Moo is large compared to 1, Ms is substantially less than 1 and 
the detached shock wave is located a relatively short distance ahead of 
the body surface (i.e., (xg + x^)/Ri 3 < < l) . The flow between the shock 
wave and the body surface is divided into two regions. Region 1 is taken 
as a domain of essentially nonviscous, non-heat-conducting, incompressible 
flow while region 2 is taken as a domain of very low speed, but compres- 
sible, viscous, and heat-conducting flow. It is .anticipated further that 
in region 2 the u and v components of velocity will be very small. The 
component of velocity w due to yaw may, of course, take on rather large 
values . 

Now it may be demonstrated with equations (l) and (2) that 8^p/8y^ 
becomes relatively independent of x along the stagnation streamline in 



Region I— Incompressible, nonviscous flow 
Region 2- Low- velocity, compressible, viscous flow 


Sketch (b) 
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the limit as the disttirhed flow extends only a short distance away from 
the body. Inasmuch as this is the type of flow of interest here, it will 
be assumed throughout this analysis that d^/Sy^ is essentially con- 
stant along the stagnation streamline between the shock and the body . 

With these assumptions, the derivative with respect to y of the 
y momentum equation yields a differential equation that becomes tractable, 
both in regions 1 and 2, when terms that vanish in the neighborhood of the 
stagnation streamline are dropped. Approximate solutions to these simpli- 
fied y momentum equations are found for the u velocity along the stag- 
nation streamline in region 1, and for the derivative of this velocity 
along the stagnation streamline in region 2 . The constants appearing in 
these solutions are determined by matching the boundary conditions at the 
shock wave and at the surface of the body, and by matching flow conditions 
at the Interface. This procedure fixes the locations of the shock wave 
and interface relative to the body . 

The energy equation is simplified in an analogous manner, and solu- 
tions valid in the neighborhood of the stagnation streamline are found 
for regions 1 and 2. The rate of heat transfer per unit area to the 
stagnation region of the body follows from the solution to the energy 
equation for region 2. 

Let us see how these thoughts apply in the case of a two-dimensional 
stagnation region. 


Heat Transfer to a Cylindrical Stagnation Region 


Zero yaw.- This problem has been treated for incompressible flow by 
Howarth (ref. 7) and more recently for the compressible flow by Cohen 
and Reshotko (ref. 8). One reason for re -investigating the matter here 
is to obtain compressible flow solutions which can be extended with rela- 
tive ease to the case of a yawed cylinder. In addition it was desired 
to obtain solutions which may be better suited to account for real gas 
effects, such as dissociation. 

To proceed, then, the stagnation streamlines are taken to lie in the 
x-z plane . The origin of the coordinate system is at the interface 
between regions 1 and 2, and the shock -wave and body-surface locations in 
this plane are Xg and x^, respectively (see sketch (a)). For the case 
of zero yaw, the z component of velocity and all derivatives with 
respect to z are, of course, identically zero. 

First a solution will be found to the steady- state y momentiom 
equation near the stagnation streamline in region 1. Since the flow is 
assumed incompressible and nonviscous in this region, equation (2) 
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simplifies to 


u^+v^ = - - ^ 
bx by P by 


(7) 


Differentiating equation (7) with respect to 


y there is obtained 


u ^ iv; + + v ^ 

bj\bxJ Sy dx \^yj ^y2 p ^^2 


(8) 


On the stagnation streamline v is identically zero and, therefore, Sv/Sx 
is also zero. In addition, the continuity equation (eq. (4)) becomes, for 
incompressible, two-dimensional flow 


|i . = 0 (9) 

OX oy 


Using this information with equation (8), one obtains 


-u 


8x^ \^xj 




( 10 ) 


Treating 8^p/8y^ as a fianction of y only, and noting that equation (lO) 
becomes a total differential equation along a line y = constant, yields a 
general solution for velocity along the stagnation streamline 


u = Ae^'' + (ll) 


where the constants A, B, and C are related by 

ItABC^ . i — (12) 

"ay" 


Note that the constants may be real or imaginary, depending on the bound- 
ary conditions . 
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Now it is anticipated that the velocity u will very nearly vanish 
at the interface x = 0 (i.e., in the sense that Uq/us < < 1); hence B 
will be approximately -A, and the corresponding approximate solution for 
velocity is^ 


u = 2A sinh Cx (13) 


To the same order of approximation, the second derivative of velocity at 
the interface, Uq'S also vanishes. The product 2AC is just the velocity 
derivative at the interface and can be evaluated from equations (lO) 
and (13)^ thus 


2AC = Uo’ = ± /- i ^ (1^) 

V P Sy 


Note that the negative root correctly describes the flow in the coordi- 
nate system of sketch (b) , since velocity decreases with increasing x. 

Consider next the steady-state y momentum equation near the stagna- 
tion streamline in region 2. In this domain viscous terms must, of course, 
be retained and thus the derivative of equation (2) with respect to y 
yields 



dp dv dv du 

u — — + p + 

dy dx dx dy 


•d^v dp dv 

dy^ ^y dy 



. 2 il 

dy^ 3 dy^ 




+ 





(15) 


Now close to the surface of the body the left-hand side of this expres- 
sion is negligible and the right-hand side simplifies so that the equa- 
tion may be written (see Appendix A) 


A f ^ 

dx \ dx^/ dy^ 


(16) 


^ Tn the limit of zero boundary-layer thickness , this solution is 
exactly the one to which equation (ll) reduces. 
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Along the stagnation streamline this eq^uation integrates to 




X + 


D 


(IT) 


The constant D is zero since = Ug" = 0 at the interface (x = O) . 

Near the surface of the body, eqixation (17) can be integrated to obtain 




ho^o * 


(18) 


In order to satisfy 



the boundary condition at the body surface 
it follows from equations (18) and (l^) that 



(19) 


Now p and 8^p/8y^ can be evaluated at the shock wave since both are 
considered constant throughout region 1. In Appendix B it is demonstrated 
that for two-dimensional flow 


' ^ uJ 

^ ( 7 ^ - i)r/ 


( 20 ) 


where Rg is the radius of curvature of the shock wave in the stagnation 
region. Substituting eqiiation ( 20 ) in equation (19) we obtain 


^b _ /Ts ~ ^s\ ^ 2 

Rb ” V 6 y Rb y 


( 21 ) 


where Re„, is the free-streaia Reynolds number based on twice the 

radius of curvatiire of the body at the stagnation point. Note also that 
the effective value of 7, the ratio of specific heats, at the shock wave 
is allowed to vary from the free-stream value. In this way, changes in 
internal molecular energy which are manifest at the high temperatvires 
encountered in hypersonic flight can be considered. 
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There remains the problem of solving the energy equation. In 
region 1, the energy eqxjatlon is simplified by neglecting all the viscous 
and heat-conduction terms . Then, for the two-dimensional problem con- 
sidered here, equation (5) reduces to 


u 


Sx 


r St _ 


0 


( 22 ) 


for which the solution is 

Tt 

CpdT = ^ (23) 

It can be seen from equation (23) that the interface temperature Tq is 
approximately the stagnation temperature T^, since the velocity at the 
interface nearly vanishes. The stagnation tempera tiire is, of coiarse , 
given by the integral equation 




7ooRW=o^ 

2 


(24) 


where for very high velocity flow the lower limit of the integral will be 
neglected. 

Next consider the energy equation in region 2. Proceeding in a man- 
ner analogous to that used in studying the y momentum equation in this 
region, we neglect the terms with the factors u, v, du/Sx, and 3v/Sy. 
Thus equation (5) becomes simply the heat-conduction equation 


bx 





(25) 


The coefficient of thermal conductivity, k, is considered a known function 
of temperature (pressxare is essentially constant) . Thus a new function of 
temperature, t), may be defined such that 


n = 



k dT 


Then eqiiation (25) may be expressed in terms of the function q 

0 


(26) 


( 27 ) 
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Inasmuch as the body boundary is cylindrical, it is convenient to use the 
general solution to equation (27) in terms of the polar coordinates (f, 0). 
Thus 


T]=A + Blnf + 


n=i 


^ ) cos n0 + 


EnF*^ + -^^sin n0 


(28) 


The origin of the coordinate system is now taken as the center of curva- 
ture of the body, and 0 as the acute angle between the radius vector f 
and the stagnation strea^ine . If a surface temperature is assumed inde- 
pendent of the angle 0, the solution on the stagnation streamline (0= O) 
reduces to 


T| = Ti, + B In — + 
•b Rb 


Cn? 




n=i 



(29) 


Letting — = 1 + e, where e is very small compared to unity, and 
Eb 

expanding equation (29) in a series of ascending powers of e, we obtain 


■^ = % + 


G e 



+ 0 (e3) 


(30) 


where G is the constant B + 


2nR^^Cn 


It is indicated by*this 


n=i 

equation that r\ varies essentially linearly with since e^/2 is 
negligible compared to e and terms of higher order in e should be 
very small indeed.^ Since e = (x-j^ - x) / r-^ < < 1^ equation (30) can be 
written 


Tl = Tlo - 




X 

^b 


(31) 


^The dependence of surface temperature on G should be small in the 
stagnation region. 

^It should be pointed out that this argument hinges implicitly on 
the assumption that T[ is a weak function of 0 near the stagnation 
streamline . 
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According to this expression, the rate of heat transfer per unit area to 
the stagnation region of the body is 




^o - % 



(32) 


The stagnation-line coordinate is substituted from eq^uatlon (21 ) , 

and the rate of heat transfer becomes 




1/4 


L8(7s - 1) J R; 


1/2 


Re, 


1/2 


Rb 



( 33 ) 


A Nusselt number is defined for interface temperature conditions 
using a characteristic length equal to twice the radius of curvature of 
the body and a temperature potential of (Tq - T^); thus 


Nu 


ko(To - Tb) 


(3i^) 


or, substituting from equation 


(33) into (3^) 


Nu 



o RP 1/2 

dpo Rb j °° 

Rs/ ko(To - Tb) 



( 35 ) 


For a relatively cool body in hypersonic flight, it is possible to dis- 
regard the lower limit of the integral and the value of body tempera- 
ture Tb compared to the interface tempera t\ire Tq. 

Note that the solutions given by equations ( 33 ) add ( 35 ) can be used 
for the case where viscosity, thermal conductivity, and specific heat are 
arbitrary functions of temperature. For instance, these functions can be 
calculated to Include the effects of vibrational and dissociational molec- 
ular energy if the extent to which these energy modes are excited is known 
throughout the flow. It is sufficient for the most part in this paper to 
consider the specific heat a constant and the viscosity and thermal 
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conductivity as proportional to the nth power of temperature. In this 


case 


( T T \ 

note —2. Rs 1 —2 > > 1 ) 

Tt J 


Too V2Cp 


7^R 


KT 


(36) 


Noting that 


1 


.Tr 


T, 


k dT = 


ko(To - Tb) To - Tb 


T> 


o , s n 

dT = ^ 
■*• 0 / n + 1 


l-(Tb/T„ 


n +1 


1 - Tb/To 


(37) 


and that 


}^oJ 






V7^R. 


M, 


-n 


( 38 ) 


it is seen that the expression for Nusselt number (eq. (35)) becomes 


Nu = 


/ xn+i 


n + 1 V^s - T y Vhsy \7c 


mJ' 1- (Tb/To) 


and the rate of heat transfer per imit area to the stagnation region of 
the body is, in terms of free-stream conditions, 




3 

8(7s - 1 )- 


1/4 



7^ 

.2Cb 


S+i 


1/2 „ n+2 ^ 

R^oo Mqo 

n+1 



(1^0) 


These considerations coii 5 )lete the zero-yaw ajialysis. However, before 
imdertaking the study of effects of yaw on heat transfer it is appropriate 
to make a few remarks. By way of introduction it should be pointed out 
that a procedure quite analogous to that just described for treating a 
cylindrical stagnation region can be employed to treat an axially symmetric 
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stagnation region. The results of such an analysis are presented in 
Appendix C • There is then the general question of the legitimacy of the 
several assumptions xmderlying the present treatment of stagnation-point 
flows. In order to shed some light on this matter it was undertaken in 
Appendix D to examine the solutions obtained to see whether they are con- 
sistent with these assumptions and with pertinent results obtained by 
others. In this regard it is shown that the presumption of a constant 
second derivative of pressure normal to the stagnation streamline yields 
solutions for the distance between shock wave and body which are quite 
close to observed values. Next^ it is demonstrated that, as assumed, the 
velocity u is negligibly small throughout region 2 \mder continuum flow 
conditions. Then it is shown that the largest of the viscous dissipation 
terms neglected in the energy equation for region 2 is indeed small com- 
pared to the heat-conduction terms. It is found too that the analysis 
predicts an amount of heat convected into region 2 which is the proper 
order of magnitude to account for the heat transferred to the body. 
Finally, it is shown that under comparable conditions equation (35) of 
this paper predicts essentially the same heat transfer as references 7 
and 8. 

In view of these results it would seem that the simplified analysis 
presented here for stagnation-region flows is, while on the one hand cer- 
tainly approximate , on the other hand quite capable of predicting useful 
information. Accordingly, we proceed to the study of effects of yaw on 
heat transfer. 

Yaw.- In this case the x direction is normal to and the z direc- 
tion IF“parallel to the stagnation line of the body (see plan view, 
sketch (c)). Then the z component of velocity has a finite value, but 
all z derivatives are again zero. 



Region /— incompressible, nonviscous flow 
Region 2^ Compressible, viscous flow 

Sketch (c) 
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The y momentum equation in region 1, differentiated with respect 
to y, takes the same form as equation (lO) on the stagnation streamline. 
Thus the velocity u is again given by the solution 

u = - i /- i sinh Cx (kl) 

C P ay^ 

The z momentum equation for the stagnation streamline in region 1 
becomes, on dropping the negligible terms from equation (3) , 


u 



(h2) 


which has the solution 


w^ = 7 RT M^^sin^A 

OO OO 00 


m 


since the transverse component of velocity is unchanged on passing through 
the shock wave . 

The energy equation for the stagnation streamline in region 1 reduces 
to a form similar to equation (22) 


u 


Sx 


+ w 


Sx 


P St _ 


0 


m 


which has the solution 



u^ + 7ooRTA^sin^A 
2 


(i^5) 


where again the stagnation temperature T-j. is given by equation (2k) . 
At the interface where the velocity u is negligible, the temperature 
Tq(A) is given by the solution to 



2 


( 46 ) 
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which, for a constant heat capacity, Cp, is 


To(A) 


The differentiated y momentum equation for region 2 takes on the 
same form on the stagnation streamline as equation (l6) . Hence, the 
solution is 




M„,^cos^ 




n ^ 2£f. 

bx by^ 2 




1 


(M) 


and the body stagnation point coordinate is 


[xi3(A)]^ = 


2Uo(A) 

7'^ dy2 


(49) 


Now, however, the second derivative of pressure is a fimction of the 
angle of yaw (see Appendix B), 



(7s - l)Rs^ 


(50) 


so the stagnation-point coordinate is given by 


^bO^) ^ ( 7 ^ - ly^^ / ^Pq(A) 2 

Rb V 6 y V Poo Rb/ 


(51) 


In region 2, the solutions to the z momentum equation and the energy 
equation are considered simultaneously. The z momentum equation simpli- 
fies (to the order of this analysis) , in the region of the stagnation 
streamline , to 



( 52 ) 
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Similarly, the energy equation near the stagnation streamline in region 2 
may he written (note that dw/dy is zero by symmetry) 





( 53 ) 


In order to facilitate the solution of equation (52), it is helpful 
to observe that the yawed boimdary layer, identified with the w compo- 
nent of velocity, resembles the boundary layer on a flat plate. It might 
be anticipated then that, just as in the case of the flat plate, the 
variation of w with x is relatively Insensitive to variations of p 
with X. In this event equation (52) has the approximate form 




( 5 ^^) 


The solution is taken in polar coordinates in order to conveniently fit 
the boundary condition that w is Identically zero at the body surface. 
Then following the same arguments used in deriving equations (29) and (30), 
one obtains on the stagnation streamline 


w = B In 


00 




n=i 



( 55 ) 


where again e = (f/R] 3 ) - 1 < < 1. If second order and higher terms in e 
are neglected, the z component of velocity on the stagnation streamline 
becomes, in terms of x/x-j^. 


W = Wq 




( 56 ) 


whence 


dw _ ^ 

dx x^ 


( 57 ) 
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If this resiilt is substituted for the last term in egviation ( 53 ) ^ 
the energy eqmtion becomes 





= 0 


m 


A solution for equation (58) which satisfies symmetry conditions on the 
stagnation streaonline and also the boundary conditions that q and n are 
constant along the surface of the body is 


Ti(X) = Tib + B In — + 
Rb 


Cnf’^ 


n=i 




where p. is a mean value of n in the stagnation region. If equa- 
tion (59) is expanded in terms of e , t] takes the following form on the 
stagnation streamline 

nU) -ni,* j(c + ( 60 ) 


The constant J is evaliiated by letting t] be tIq when 6 is Cq = Xb/Rt 
and is given by the relation 


J = 




2eo^ 


(1 + 6 q ) + . . . 


(61) 


The rate of heat transfer per mit area to the stagnation region of the 
body at angle of yaw X is, from eqxxation (60), 


8n(x) 

8x ib 


8ti(X) 
Rb 


8e 




26n^ 


(62) 


Substituting equations (26) and (61) into this expression and neglecting 
ter*ms of the order Gq compared to 1, one obtains 


-nb'(A) 


1 

Xb(A) 



k dT + 




( 63 ) 
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Multiplying by 


gBh 


and substituting from equation (51) yields 


gnh'(^)Rh 

^o(To-Tb) 



Vm-o(A) ^s/ 


1/2 1/2 

Repp cos A 

ko(To - %) 



(64) 


For a constant heat capacity it follows from equations (23)^ (43) > 
and (47) that 


To(?^) 

To 


To 




cos^A 




2CpSin^A 


J 


(65) 


If, in addition, the thermal conductivity is proportional to the nth 
power of temperature, then 


^0(^0 " ) 


pTo(A) 

/ k dT = 

U m 


2n+2^ 
COS A 


(n + 1)(1 - Tb/To) L 


1 - 


TqCOS^A 


n+1 


and 


( 66 ) 


dWo^ 


2ko(To - Tb) 


p[ \Prsin^ 
.^o A - Tb/To 


T. 2 n-> , 2, 

b Pr cos A sin A 

1 - Tb/To 


(67) 
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Thus equation (6^4-) becomes 




1 




ko(To - %) n + 1 Vs • V \W M^^(i _ t^/t^) 


{ 


cos^A 


1 - 


Tb 


TqCOS^A 


n+i 




( 68 ) 


The ratio of equation ( 68 ) to equation (39) is the ratio of* the rate 
of heat transfer to the stagnation region of a yawed body to the rate of 
heat transfer to the stagnation region of the same body at zero yaw. This 
ratio is 


X'M f ,, 

n^'(o) l-(Tb/To) i 

( 69 ) 

An analogous expression can be obtained for the ratio of Nusselt 
numbers , thus , 


1 - 


TqCOS A, 


+ (n +l)Pr 




sin^A 


I 


Nu(A) ^ %'(A) Tq - 

Nu(o) %'(o) Tj. - Tb 


(70) 


where from eq.uation ( 63 ) the recovery temperature, Tj-j is the solution to 


k dT = - ( 71 ) 


However, it should be noted that the ass\jmptions used in the analysis 
tend to be violated when the body temperature approaches recovery condi- 
tions. Therefore it should not be expected that equation ( 7 I) will yield 
accurate values for recovery temperature. 
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There remains, of coiirse, the problem of determining {I. It is suf- 
ficient for the pxrrposes of this report to take p as s imp ly the arith- 
metic average between iio(^) and that is, ji= [hq(A) + ii^]/ 2 .^ In this 
event equation (69) can be written 


1b'(o) 


n+a/ 2 , 

cos A 


1 - (Tb/To) 


n+i 


cos^A 



+ 


n+ 1 
2 


Pr 



• 2, 

sin A 


} 


(72) 


which in the case of a relatively cool surface (i.e., T-j^/t^cos^A < < 1) 
becomes 


%'(A) 

nb'(o) 


n+ 1/2 

cos 


A 



n+ 1 
2 


Pr 


sin^A 


(73) 


These considerations complete the theoretical analysis. Attention 
is turned next to the experimental investigation. 


EXPERIMENT 
Test Apparatus 


In order to investigate hypersonic heat-transfer phenomena, it was 
undertaken to develop an apparatus capable of producing a hot hypersonic 
air stream. In addition to the air stream being hot (i.e., characterized 
by stagnation temperatures of the order of thousands of degrees Rankine) , 


4 

Act\ially this procedure might better be considered the first step in 
an iteration method where fl is recalculated on the basis of the preced- 
ing calculation of T as a function of x. This refinement is not con- 
sidered warranted here where only the gross effects of yaw for angles of 
yaw well less than 90 ^ are of principal interest. As the angle of yaw 
approaches 90°, the analysis as a whole tends to break down due to the 
violation of the several assumptions predicated on the flow being hyper- 
sonic normal to the axis of the cylinder. 
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it was required that the stream he of sufficient duration to facilitate 
measiirements of heat transfer. The apparatus which finally evolved is 
termed a hypersonic gun tunnel simply because it combines a gun with a 
wind tunnel. Inasmuch as this gun tunnel has not been described in gen- 
erally available literature, a brief description of its operation, cali- 
bration, and performance will be given here. 

A schematic diagram of the equipment is shown in figure 1. A smooth- 
bore gun is used essentially as a one-stroke pimp to create a reservoir 
of high-press\ire, high-temperature air. A charge of fast burning powder 
is biirned at the breech to accelerate the light-weight piston to a high 
velocity in the barrel. The accelerating piston produces strong shock 
waves which reflect a nimiber of times from the blanked end of the barrel 
and the piston. These shock waves compress the air in the gun barrel in 
a highly nonisentropic manner, transforming much of the energy available 
in the powder charge to the form of heat in the air. Accordingly, a res- 
ervoir of very high-temperature air should be created.^ As soon as 
equilibrium pressures have been reached in the reservoir, the valve clos- 
ing the entrance to the nozzle is opened and the high-ten^erature , high- 
pressure air expands through the nozzle into the test section where models 
are located. 

Figure 2 is a photograph of the present model of the hypersonic gun 
tunnel. This unit consists of a 20-mm smooth-bore barrel 5 feet long 
with a breech chambered for the standard 60/20 cartridge. The barrel is 
initially charged with dry air at 10 atmospheres pressure and room tem- 
perature.® A 20-mm dieuneter piston constructed of nylon and weighing 
3-3/^ grams is used, and the cartridge is loaded with 27 grams of IMR 
No. ^227 rifle powder. The nozzle valve is released 80 milliseconds after 
the powder charge is ignited. In this time the reservoir conditions, as 
indicated by pressure measxorements , have approached steady values. The 
flow exhausts through the nozzle into a vacuum tank. This tank, seen in 
figure 2, has a volume of 10 cubic feet and is pimrped down to a pressure 
less than 0.0007 atmosphere before each test. In this manner sufficient 
compression ratio (up to ^5X10^) is provided across the nozzle to maintain 
hypersonic flow in the test section during the course of the run. 

The nozzle through which the flow expands is fabricated of stainless 
steel and has a simple conical contoiir of 9° total angle and a 1-inch 

^In Appendix E an estimate is made of the final temperatures that 
may be obtained from this nonisentropic compression. It is indicated 
that for a compression ratio of about 32, the ratio of final to initial 
temperature would be between 4 and 7. The corresponding isentropic com- 
pression process would yield a final to initial temperature ratio of only 
about 2. 

^Higher final air temperatiores sho\ild, of course, be obtained if the 
initial air temperatiires are increased. This possibility is presently 
being investigated. 
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exit diameter (see fig- 3) • While the flow is thus expanding conically, 
the rate of expansion is slow in the region of the models and for most 
practical purposes the flow may be considered parallel. The throat of 
the nozzle is a cylindrical section 0.020-inch in diameter and 0.050-inch 
long. This section is machined as a separate insert and can be replaced 
when the nozzle throat becomes seriously eroded by the heated air- 
Because of the length of the throat section, several runs can usually be 
made with one insert before erosion becomes serious in the sense that it 
increases the effective sonic throat diameter. 


Calibration of Flow 


Properties of the hypersonic stream produced by the gun tunnel were 
obtained from the following four measurements: reservoir pressure, pitot 

pressure, static pressure, and stream velocity. The reservoir pressure 
was measured with the gage shown in figure ^(a) . The piston of the gage 
was housed in an insert which connected to an orifice leading to the 
reservoir section of the barrel (figs. 1 and 2). The 0-ring seal near 
the bottom of the piston prevented leakage from the reservoir to the 
atmosphere, and the entire unit was retained by a yoke which was bolted 
to the side of the reservoir. The piston transmitted the pressure force 
from the reservoir to strain-gage elements which were electrically con- 
nected to a conventional bridge circuit, amplifier, and recording oscil- 
loscope. The time required for the reservoir pressure to reach equilib- 
rium was deduced from measurements made with a low-sensitivity gage 
element. A typical record of the initial pressure pulses is shown in 
figure 5(a) , and it can be seen that these pulses are damped in well less 
than 0.06 second. The measured peak values of the first pulse or two are 
probably lower than the actual pressure maxima because of the inertia of 
the gage.^ 

The reservoir pressures used for calibration were measured with a 
high-sensitivity-gage element which was mechanically stopped at 400 atmos- 
pheres so that the peak pressure would not damage the gage. The average 
pressiire for all runs is shown in figure 5(b) together with the standard 
deviation from the average. This deviation was within ±10 percent. The 
timing pulse shown in the figure corresponds to the signal superimposed 
on the oscilloscope record when the nozzle valve was opened. This pulse 
was taken as the zero time reference for the beginning of flow through 
the nozzle . 

The pitot pressure was measured with a conventional, strut-mounted 
probe constructed from a 0.050-inch O.D. hypodermic needle (fig. ^(b)) 
with an outer to inner diameter ratio (O.D./l.D.) of 1.47. The probe led 

It may be of interest that the frequency of these measiired piilses 
is nearly identical to that predicted by the method of characteristics. 
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to a variable-capacitance diaphragm-type pressiare cell (fig. ^(b)) con- 
nected to a bridge circuit, amplifier, and a recording oscilloscope. 

The response time of this system appeared to be instantaneous for prac- 
tical purposes. The average value of the ratio of the pitot pressure 
(measured on the nozzle center line) to reservoir pressiire is shown in 
figure 5(c). This ratio is approximately constant at 6.6x10”'^ during 
the first 0.3 second but then, for some unexplained reason, jumps to a 
higher value and holds approximately constant imtil about 0.6 second. 

For the first 0.3 second, standard deviation limits shown in figure 5(c) 
are within 4 percent of the average. 

The static pressure was also measured with a conventional type probe 
consisting of a 13 ° included angle cone-cylinder combination of 0.050-inch 
O.D. mounted in the same manner as the pitot probe. The pressure orifices 
were located at the test station and, to determine if the measured values 
of static press\ire were grossly influenced by Reynolds mmiber effects or 
feedback, a series of probes were tested. These probes had forebodies 
and afterbodies that ranged from 30 to 6o diameters in length as measured 
from the pressure orifices to the tip of a probe and to the leading edge 
of the strut, respectively (fig. ^(c)). Each probe was connected to a 
capacitance pressure gage which differed in only one respect from that used 
for pitot pressure measurements. This difference was that a thinner dia- 
phragm element was used in order to obtain the necessary instrument sensi- 
tivity at the lower static pressures. The pressures measured on the nozzle 
center line with the several probes were essentially in agreement, and the 
value of measured static to reservoir pressure for all runs is shown in 
figure 5(d). The standard deviation limits for the first 0.3 second are 
within 8 percent of the average . 

Now it seems reasonable to expect that the stream static pressure 
actually varies like the pitot pressure during this interval of time; 
that is to say, changes in p/H with time (fig. 5(b.)) are due primarily 
to the finite response time of the static pressure measuring system. 

With this thought in mind the measiired data on p/H have been reduced 
in the manner described in Appendix F . The corrected value of* P/H is 
found to be 5.32x10"® for the first 0-3 second of rtm. The corresponding 
value of p/pt diuring this interval is 8 . 06 xl 0 “®. 

Pitot and static pressures were also measured l/8-inch off the nozzle 
center line, and they were found to agree with the center-line measiire- 
ments to within the standard deviation. 


®The flow does not actually choke to subsonic speeds until about 
1.0 second, but after 0.6 second it is rather erratic. Accordingly, the 
flow was calibrated only over the interval from 0 to 0.6 second. It was 
fovind that the jump in pitot pressure occvirring between 0.3 and 0.4 sec- 
ond is accompanied by unpredictable changes in static pressure, in shock- 
wave profiles about models, and in heat-transfer measurements. 
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The stream velocity was obtained by determining the rate at which a 
disturbance passes through the test section. This disturbance was pro- 
duced by the spark discharge of a 0-5 mfd capacitance at 10 kilovolts 
across electrodes mounted in the stream ahead of the test region (see 
fig. 4(d)). The position of the disturbance in the test region was 
recorded photographically with a schlieren apparatus using a spark dis- 
charge light so\nrce. The velocity of the stream was then taken as the 
distance between the tips of the spark electrodes and the center of the 
disturbance divided by the time interval between the disturbance spark 
discharge and the schlieren exposure.® Some typical schlieren photographs 
of different disturbances are shown in figure 6(a) . (At the bottom of 
the test region can be seen a blunt test probe and its associated shock- 
wave system.) Measurements of velocity were made for separate runs over 
the range of flow time from 0.03 "to 0.6o second. The flow time was again 
measiored from the opening of the nozzle valve. At each flow time, a 
series of runs were made in which the time interval between the disturb- 
ance discharge and the schlieren exposiire was varied. The data, shown 
plotted in figure 6(b) , appear to be relatively independent of this time 
interval, and it is therefore concluded that any consistent error exist- 
ing in the velocity measurements is small. 


The reservoir pressure was not used directly to determine stream 
properties, since it was considered unjustifiable to assume that the 
expansion of air from reservoir to test region satisfied the usual isen- 
tropic relations. However, the reservoir pressure was used to nondimen- 
sionalize the meas\ired static and pitot pressures, inasmuch as the ratios 
p/h and P-j^/h were more consistent from r\m to run than the absolute 
pressures . Pertinent stream properties were then determined in the fol- 
lowing manner. At hypersonic Mach numbers, the ratio of pitot pressure 
p^ to dynamic pressure is approximately independent of Mach number, 

namely. 


7 7+1 



The function P( 7 ) is very nearly constant, varying only from 0.920 to 
0.977 a-s 7 varies from 1.4 to 1.1. Then, the unknown stream parameters 
are, in terms of the measured properties 


density, 


static temperature, T„ = = (FU„^^)p^p^ ^ 

Mach number , M^ = 


(75) 


®The electrodes were of slender ( 15 ° wedge 0.075 inch thick at root) 
cross section in the streamwise direction, so that while they might alter 
stream Mach number, they were unlikely to appreciably influence stream 
velocity at the hypersonic speeds of these tests. 


CONFIDENTIAL 


MCA EM A55E02 


COHFIDENTIAL 


29 


and the stagnation temperature is the solution to 

/.Tt 2 

ip = 2 

The measured velocity of the stream is nominally 5000 ft/sec. 

Then, if ideal gas values are used for heat capacity in equation (t 6) , 
the stagnation temperature is 2200° R. If vibrational energy modes are 
fully excited in any part of the flow, the correction for this caloric 
imperfection might decrease the stagnation temperature to 2100° R (see 
refs. 9 and 10). This correction is within the deviation in velocity 
data, however, and therefore was not applied. 

The ratio of static to pitot pressure, 8.o6xlO“®, yields a Mach num- 
ber of 9.8 (eq.. (75)) and the initial values of the pitot and static pres- 
siire are 0.22^ and O.OOI8 atmosphere, respectively. These pressures, 
and therefore the density also, fall off approximately exponentially dur- 
ing the flow. Accordingly, the nominal free-stream Reynolds niunber also 
decreases with time and is given approximately by the relation. Re = 1.26 
g-t/0.6 million per foot. 

The stream properties given above are all reproducible within 10 per- 
cent. It will be noted that the stream velocity appears somewhat lower 
than the nominal value at the beginning of the flow (see fig. 6 (b)). This 
result may be due to the initial cooling effect of the sonic throat region 
and the reservoir. In any case, however, this decrease in velocity is 
within the standard deviation of the data, and the nominal value of veloc- 
ity will accordingly be taken as 5000 ft/sec for the entire interval of 
flow. 


In addition to the above calibration, the piirity of the air stream 
was tested by analyzing gas samples with a mass spectrometer and samples 
of condensable phase impiirities by X-ray diffraction and emission spectro- 
graph methods. The gas samples were collected with a pitot tube leading 
to an evacuated glass flask. At a predetermined time of flow, the copper 
tubing leading to the flask was sealed with a hydraulic pinch. The con- 
densable phase impurities were simply collected by impact on a spectro- 
graphic grade copper rod and time scale effects were differentiated by 
using a movable rod within a shield. 

It was found that the gas samples were essentially normal air except 
for about 1- percent carbon dioxide and a trace of nitrous oxide. As 
expected, the percentage of these impurities increased as the sampling 

^°Velocities up to about 80OO ft/sec have been obtained using larger 
powder charges. This result implies stagnation tenperatures in excess of 
5000° R. Damage to the gun barrel was, however, excessive to the point 
of being intolerable under these circijmstances . 
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time interval was lengthened. The condensable phase ingjurity was found 
to be iron oxide, undoubtedly originating from observed areas of erosion 
in the reservoir and the nozzle throat. Samples taken at the beginning 
of flow showed very light film deposits of impurities, and for this reason 
it is felt that the initial heat-transfer rates measiired at the beginning 
of flow shoiild most nearly represent the correct values for nonnal air.^^ 
Therefore, only these initial data are interpreted in the following por- 
tions of this paper. This procedure also serves essentially to eliminate 
the problem of evaluating heat-loss corrections for the heat-transfer 
models since these corrections must vanish at zero time when thermal gra- 
dients along a model are zero. 


Heat-Transfer Models and Reduction of Data 


The models used in the heat-transfer experiments were tested at the 
same station at which the pitot and static pressure data were taken. 

These models were butt-welded iron-constantan thermocouple cylinders 
(fig. 7(a)) transverse to the flow at various angles of yaw. Wires of 
0.003- , 0.010-, 0.012-, 0.020- , and 0.0^0-inch diametef were tested at 
zero yaw. The 0.003- inch-diameter wire was tested at 22.5°, and 70° 

angle of yaw and the 0.020-inch diameter wire was tested at ^5° yaw (shown 
installed in the test region, fig. 7(b)).^^ The thermocouple junctions 
were examined microscopically and were rejected if the weld showed visible 
imperfection. 

The diffusivity of the thermocouple material was sufficiently large 
to permit the assumption that ten5)erat'ure gradients across the thermo- 
couple junction were negligible (i.e., the time of diffusion to within 
95 percent of constant temperature is very much less than the time con- 
stant of response of the thermocouple to heat transfer from the air 
stream) . Thus the average rate of heat transfer per unit area to the 
models was proportional to the rate of change of temperature indicated 
by the thermocouple junction. Under these conditions the average heat- 
transfer rate satisfies the relation 


- , I^bCbPb 

\ ~ 2 dt 


( 77 ) 


^^The possibility of these impurities causing errors in the heat- 
transfer data is considered fiirther in Appendix F. 

^^The ratio of length to diameter of all test cylinders was 25 or 
greater. In the case of the yawed cylinders this ratio was always in 
excess of 50. Accordingly, it is felt that end effects on measured heat- 
transfer rates are negligible, with the possible exception of the data 
for 70° yaw. 
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where is the temperature of the thermocouple with radius speci- 
fic heat C-t), and density p-|^ (note is essentially the same for 

iron and constantan) . It is convenient in the analysis of data for the 
cylinders at zero yaw to employ an average heat-transfer coefficient h 
defined by the equation 



(78) 


It follows that the average Nusselt n\amber, based on stagnation tempera- 
ture conditions and a characteristic length eqijal to the diameter of a 
cylinder, is given by the expression 


= = _!b!SbPk_£^ 

l^t kt(Tt - Tb) 


(79) 


Figure 8(a) shows the output from the 0.003-inch-diameter transverse 
thermocouple for a typical zero-yaw run. The rate of temperature rise 
drops below expected values after about 0.05 second. At about 0.3 second 
the record indicates some perturbation in the flow as did the static and 
pitot pressure measurements, and observed shock-wave profiles. The initial 
temperature rise, which is of principal interest here, repeated to within 
±10 percent from r\m to run. In order to measure this initial rise more 
accurately, the oscilloscope sweep speeds were adjusted to give about 45^ 
slope of the temperature -time curve. A typical temperature -rise curve on 
the expanded time scale is shown in figure 8(b) . 


RESULTS AND DISCUSSION 


The theory of this report treats, of course, only the heat transfer 
to the stagnation region of a body in hypersonic flow. Nevertheless, for 
a relatively blunt body, such as a transverse cylinder, it seems reason- 
able to expect that in hypersonic flow, just as in lower speed flow, the 
average rate of heat transfer will vary in proportion to the stagnation- 
region heat-transfer rate. With this thought in mind, it is noted that 
according to equation (39) ^*0^ ^ cylindrical stagnation region 

(A = 0, n = 1/2, Tb/To < < 1) 
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If R^/Rg is considered constant (presimably at some value near l ) , 
then it is easily demonstrated that 

Nu ~ /S ( 81 ) 

neglecting small differences between m-q and |j.o . This is essentially 
the result deduced experimentally by Kovaznay (ref. 11) and by others 
for the average Nusselt numbers of cylinders transverse to subsonic and 
supersonic flows. The experimental data obtained in the hypersonic gun 
tunnel should provide some indication of whether or not the result holds 
in the same sense for hypersonic flows. These data for transverse cylin- 
ders are presented in figure 9,^^along with the predictions of equa- 
tion (8l) (replacing Nu with Nu) and with data obtained from hot-wire 
type experiments by Kovaznay (ref. ll) , Stine (ref. 12), and Stalder, 
Goodwin and Creager (ref . 13) • Specifically, there are shown (after the 
method of Stine) the average Nusselt nimibers as a function of M yRe^Pr. 
The agreement between the data obtained from the hypersonic gun tunnel 
and those obtained by the hot-wire method is surprisingly good consider- 
ing the wide difference in test conditions'^ (e.g., as contrasted to the 
gun-tunnel data, the data of reference 12 are for Mach numbers less than 
1.4 and stagnation ten^eratures less than 590° R) . Note, too, that the 
scatter of the data from the gun tunnel is by and large of no greater 
order of magnitude than that of the data in references 11, 12, and 13 . 
Finally, it is observed that equation (8l) is quite as consistent with 
experiment at hypersonic speeds as it is at lower speeds. Evidently, 
then, average rates of heat transfer to transverse cylinders in hyper- 
sonic flow do tend to vary in proportion to the stagnation-region heat- 
transfer rates . 

It is appropriate to inquire now about the effect of yaw on the rate 
of heat transfer to a cylinder in hypersonic flow. In this regard it is 
of interest to see first what theory predicts for the effect of yaw on 
heat transfer to a cylindrical stagnation region. Two cases will be 
treated, namely, T^/Tq = 0 and T^/Tq = 0.24, the value corresponding to 
test conditions in the hypersonic gun tunnel. If Prandtl number is taken 
equal to 3/4, and the tempera tiire exponent for the thermal conductivity 


^^Note that while the test Reynolds numbers are low, they are, accord- 
ing to estimate, within the range of continuum laminar flow. 

^"Sniile the stagnation temperatures of the hypersonic gun tiinnel were 
far greater than those of the other tests, it should be noted that they 
were well below temperatures at which dissociation might be expected to 
occur and, hence, possibly influence heat transfer. 
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and viscosity functions is taken as then according to equation (72) 

the stagnation-region heat-transfer rate is influenced by yaw in the man- 
ner shown on figure 10- For the case of T^/Tq = 0, it is indicated that 
yaw shoTild reduce the rate of heat transfer almost in proportion to cos^A 
(note this is the same reduction factor as for impact pressures) . For the 
case of Ti^/Tq = 0.2h the predicted reduction in rate of heat transfer is 
somewhat less, although still very sizable. How well this prediction 
agrees with measurements of average heat-transfer rates to yawed cylinders 
in the hypersonic gun tunnel is shown in figure 11. The average rates are 
observed to decrease with increasing yaw in much the manner predicted by 
the stagnation-region theory. For example, at 70° yaw, theory and experi- 
ment show about a 6o-percent reduction in heat-transfer rate . 

The stagnation-region theory and hypersonic gun-tunnel experiments 
have, then, provided useful information on the problem of heat transfer 
to cylinders in hypersonic flow. However, all of the considerations, 
both theoretical and experimental, up to this point have been for the 
case where the gas does not dissociate and therefore the gas properties 
vary with temperature in a somewhat ideal manner . At the very high tem- 
peratures encountered in hypersonic flight (namely, at temperatures in 
excess of about 5000° R) gas molecules may dissociate into atoms. It is 
natural then to inquire how dissociation might alter our previous con- 
clusions. Accordingly, the effects of equilibrium dissociation on the 
rate of heat transfer to a cylindrical stagnation region at zero yaw have 
been examined theoretically in Appendix G. The calculated heat -transfer 
rates are increased by dissociation at Mach numbers greater than 10 (see 
fig. 12). It is suggested, then, that equilibrium dissociation (should 
it occur) may alter oiir previous conclusions regarding the effect of yaw- 
ing a cylinder to the extent that somewhat greater reductions in heat- 
transfer rates will be achieved. This possibility derives from the fact 
that the high temperatures which bring about dissociation tend to be 
decreased by yaw.^® 


CCfl7CLUDING REMARKS 


It has been argued theoretically and demonstrated experimentally that 
yaw can have the effect of substantially reducing the rate of heat transfer 

^^The 1/2 power law is, at the very high temperatures encountered in 
hypersonic flight, a rather good approximation according to Sutherland's 
equation (Cp = const) , although actually the results of theory are rela- 
tively insensitive to n (except at very large yaw) over the range 
1/2 < n < 3 / 1 ^. 

^°It is perhaps worthy of note, too, that by the same token radiation 
heat transfer from the high temperature disturbed air to the body surface 
shoiald be alleviated by yaw. 
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to a circular cylinder in hypersonic flow. Experiments on small cylin- 
ders in a hypersonic air stream of Mach number 9.8 and stagnation 
temperatiire 2200° R tended to verify this prediction, indicating that 
the rate of heat transfer was reduced to about 4o percent of its zero- 
yaw value when the angle of yaw was 70°. These results are interpreted 
to have the practical significance that sweep may markedly reduce the 
heat transfer per unit area to the bl\mt leading edge of a wing in hyper- 
sonic flight. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
May 2, 1955 
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APPETOIX A 


SIMPLIFICATION OF THE y MCMENTUM EQUATION IN REGION 2 

The steady- state , two-dimensional y moment\jm equation (see eq. (2)), 
differentiated with respect to y, yields 



Now on the stagnation streamline the velocity v is identically zero 
and therefore all x derivatives of v are zero* Also^ all odd order 
y derivatives of functions like density p, viscosity p, pressure p, 
and velocity u vanish since, by symmetry, these functions are even. 

In addition, it is assumed that near the stagnation streamline the veloc- 
ity u is so smfl 1 1 throughout region 2 that terms with this factor may 
be neglected. With this assumption an additional useful relation can be 
deduced from the continuity equation 




u ^ + p 
dx 




+ v|£. 0 

oy 


(A2) 


Eliminating the terms with factors u, v, or ^p/By from equation (A2) 
there results, as for incompressible flow. 


^u ^ dv 
Sx 3y 


(A3) 


Note that all derivatives of the sum Su/dx + hv/hy are also zero in the 
regions where equation (A3) will hold. 

Applying the above considerations sin 5 )lifies equation (Al) to 


P 





+ pT— + — ^ 

hybx _ V^y / 


(a4) 
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Nov it will be asstimed, as is usml, that the viscous flow in the region 
of the stagnation point of a blunt body is similar to viscous flow at 
the stagnation point of a body with infinite radius of curvatiare insofar 
as the velocity derivatives are concerned (i.e., the principle effect of 
the body ciirvature is to determine the magnitude of the pressiire deriva- 
tives) . Accordingly, B®v/Sy^ and d^u/dy^ will be supposed to vanish 
in the stagnation region. Then expansion of the second member of the 
right side of equation (a 4) yields 

p 6v ^ dp 5^v d^v’ N 

6y ^ dy dy2 ^ ^ 


in which the only term retained is 2(d^ii/Sy^) (dv/dy) . Similar expansion 
of the last member of equation (a 4) gives 


(hi + + ^ + h. (h}L + + d\ ^ 

dxdy \dy dx/ dy Vdxdy dx^/ dx dxdy/ Vdy^dx dydx / 


Note that from equation (A3), d^u/dy^dx is equivalent to -(d®v/dy^) and 
will therefore be neglected. The terms retained in this equation, then, 
are (dp/dx) (d^v/dxdy) + p(d^v/dydx^) . These terms can be combined into 



Equation (a 4) thus is reduced to 


A f p A /^6vY 

dx \ dx^y dy^ dy^ dy ^ vdy/ 

The derivative dv/dy vanishes at the surface of the body, so that in 
the immediate region of the stagnation point, eqmtion (A5) takes on the 
approximate from 



^ ( d^u\ ^ _ d^ 

dx V dx^ J dy^ 


(a6) 


This expression will be taken to hold near the stagnation streamline 
throughout region 2. 
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APPENDIX B 

BOUNDA.RY VELOCITIES AND PRESSURE DERIVATIVES 


For hypersonic Mach numbers, the density ratio across an oblique 
shock wave is 



7s + 1 
^s - 1 


(Bl) 


then the pressure just downstream of the shock is 


Ps 


7s + 1 


(B2) 


where a is the acute angle between the shock wave and the normal to the 
free-stream velocity vector (see ref. 9)* I't can also be shown that the 
V component of velocity just downstreaun of the shock is 


V 


s 


2 

7s + 1 


HxjSin a cos a 


(B3) 


while the u component on the stagnation streamline is 


u 


s 


7s ■ 

7s + 1 


(bU) 


In evaluating the derivatives , consider a shock wave with radius of 
curvature Rs • Let s be the distance along this profile measured from 
the stagnation streamline and x(s) and y(s) be the equations for the 
shock-wave coordinates. Then 


while 


dv ^v dx ^ ^v dy 

ds Sx ds Sy ds 


(B5) 



p Sfp_ ^ ^ ^ dp d^x d^y 

yds J Sx^y ds ds \ds J ds^ by ds^ 


(B6) 
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In terms of the radius of curvature Rg> the differential equations for 
x(s) and y(s) are 





(B7) 




and at the stagnation streamline (ds 


= O) the following conditions hold 


ds 

dx 

di" 

d^y 

ds^ 

d^x 

ds^ 


1 

0 

0 

2^ 

Hs 








(b8) 


Then, at the stagnation streamline, eqiiations (B5) and (b 6) become 


and 



,8y^/ ds^ Rs \8x/s 


(B9) 


(BIO) 


Now by continuity and equations (B3) and (B9) 



2U„ 


(7s + l)Bs 


(Bll) 
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for two-dimensional flow. For axially symmetric flow the corresponding 
relation is 




(7s + l)Rs 


(B12) 


According to equation (B2) the first right-hand term of equation (BIO) is 


+ i)Rg 2 


(b13) 


while the next term, - -i- (— ) , is evaliiated using the x momentum 

Rs V^x/g 

equation (eq. (1)) which for the nonviscous incompressible flow region on 
the stagnation streamline reduces to 


dx 



(Bli^) 


According to equations (Bll) and (B12) , equation (B1^) becomes 


^ 2(73 " l) 

Sx (7 +1)2 Rg 


(BI 5 ) 


and 


Sp _ ^(Ts ~ Ps^oo^ 
Sx (7g + 1)^ Rs 


(bi6) 


for the two-dimensional and the axially symmetric flow cases, respectively. 
Then the corresponding second partial derivatives of pressiire are 


^ = 6 (7s - DPsUqo 

ay^ " 


( 7 , + 1 )%^ 


(BIT) 
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and 



()-s + 1) V 


(Bl8) 


Note that 7g can have values somewhat different than lA if vibra- 
tional and dissociational energies are excited at the shock wave. The 
results of this appendix are consistent if 7 g is defined by equation (b 1) 
from the ratio of densities across the shock wave. When additional energy 
modes are excited at the shock wave, this effective value of 7g is not 
exactly the ratio of specific heats . 


It can be seen that for the case of a yawed two-dimensional body, the 
same relations hold as for the body at zero yaw except that the velocity 
is replaced by the normal component of velocity, UngCos A. Thus the 
yawed two-dimensional body has a second derivative of pressure 


^ ^ _ 6(7^ - l)p^U^^cos^A 
(7s + 


(BI 9 ) 


In the above relations the radius of c\irvature of the shock wave Rg 
is yet undetermined. In the limit of infinite free-stream Mach number, 
the ratio of shock wave to body curvature, Rg/Rb^ might be expected to 
approach unity as an upper bound. On the other hand, a value of Rg/Rb 
consistent with incompressible boundary-layer solutions may be a reason- 
able lower boimd. In this regard Howarth (ref. 7 ) reports that for two- 
dimensional flow 


8v _ ^s 
8y R-b 


(B20) 


which, according to equations (B^) and (Bll) , 

Rb 7-1 


corresponds to a ratio 


(B21) 


Sibulkin (ref. l4) , using a similar analysis finds that 

8v _ 3Ug 

8y 2R^ 


(B22) 
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for axially symmetric flow. This corresponds to the ratio 


^ ^ 

Rb ~ 3(7 - 1 ) 


(B23) 
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APPENDIX C 

HEAT TRANSFER TO AN AXIALLY SYMMETRIC STAGNATION REGION 


The methods used to calculate the rate of heat transfer to a cylin- 
drical stagnation region can also be applied to the stagnation region of 
a spherical body. This analysis is parallel to that for the cylinder at 
zero yaw and thus the x axis is taken as the stagnation streamline and 
the origin of the coordinate system is placed at the interface between 
the assumed incompressible nonviscous region 1 and the viscous, low- 
velocity, compressible region 2. For the purpose of obtaining the solu- 
tions for velocity in regions 1 and 2 on the stagnation streamline, it 
is most convenient to consider the momentum and continuity equations in 
cylindrical coordinates . Because of axial symmetry, all properties are 
independent of the angular coordinate and, accordingly, the r direction 
momentum equation becomes 


P 




^r 


2 

3 Sr 




+ 




+ 


2 


r VSr 





(Cl) 


While the continuity equation is 


^ (pu) + 7 ^ (prv) = 0 


Sx 


(C2) 


In region 1 where the viscous terms are considered identically zero, 
the r momentum equation (eq. (Cl) ) reduces to 


P 


Sv 
"" Sx 


V 


Sv 

Sr 


Sr 


(C3) 


Differentiating equation (C3) with respect to r and dropping terms with 
factors V and Sv/Sx, which vanish on the stagnation streamline, gives 


/SvY , ^ S\ _ 1 s^p 

\Sr/ SxSr P Br^ 


(C4) 
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Now the continuity equation (eq.. (C2)) expands to 

p ^ + PI + p ^ = 0 (C5) 

bx r br 


however, on the axis of symmetry, neglecting terms higher than second 
order in r , 


^ = I (C6) 

dr r 

Thus, for incompressible flow, the continuity eq'uation on the stagnation 
streamline reduces to 


du dv 
dx ^ ^ dr 


(C7) 


Substituting equation (C7) into equation (C4) yields 


i. = - i. ^ P (c8) 

4 \bxj ~ 2 dx2 P dr^ 


which, upon differentiating with respect to x, and assimiing 
1 b^T) 

— = constant, becomes 

P dr^ 


u 


d®- 


u 


dx"" 


= 0 


(C9) 


For nonzero values of the velocity u, this differential equation has as 
a solution 


u 



+ Uq'x + Uq 


(CIO) 
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The value of velocity at the interface, Uq, is again considered very 
small. Thus, from eq,uation (c8) , the first derivative of velocity at the 
interface is, approximately. 


u^ ' = 



(Cll) 


As can be seen from the solution for velocity^ equation (CIO) ^ the second 
derivative of velocity is constant. Therefore the second derivative of 
velocity may be evaluated from equation (c8) using conditions just behind 
the shock wave^ thus, 




,2 

^s ^ 2 

2ug ^sPs 



(C12) 


Substituting for the values of velocity, velocity derivative, and 
second pressure derivative behind the shock vave (see Appendix B) yields 


and 


Ur 


8(3 - 27 ^) 

?'s^ - 1 Rs^ 


(C13) 


Ur 


- 1) Uqo 

7s + 1 


(Cl4) 


Now in region 2, the viscous terms are retained in eqioation (Cl) • 
Following the procedure used in studying two-dimensional flow (see Appen- 
dix a), the r momentum equation, differentiated with respect to r, is 
simplified to 


8x \2 8x^ / 8r^ 


which integrates to 


(C15) 


2 dr^ 


(Cl6) 
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and, as in the case of the two-dimensional flow, 


!£ ^ 

2 dx 


2 

o p X 
2 


+ Ax + B 


(C17) 


The constants A and B are again determined by matching the first and 
second derivatives of velocity at the interface. Thus 


A = 


^"o^o" 

2 


2 


(Cl8) 


At the body du/dx vanishes, and solving for the coordinate x-j^ from 
equation (Cl?) resiilts in 


^o^o" 





48 p Uq' 


8r^ Po^o" 


(C19) 


In Appendix B it is shown that 



8 Pcoy 

7s + 1 Rs^ 


(C20) 


Thus from equations (C13) , (Cl4), and (C20), it can be shown that 


^ ^o* ^ (7 - 1)^ n/2(7 - 1) ^ ^ 

^ 8r^ PqUo"^ " (3 - 27)2 ^^o ^^b 


(C21) 


which is large compared to unity for any reasonably large value of 
Reynolds nxamber (of the order of himdreds or greater). Therefore, if 
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quantities of the order of unity are neglected in equation (C19) ^ the 
stagnation point coordinate reduces to 





(C22) 


which is identical in form to the relation for body surface coordinate 
in the two-dimensional flow (eq. (19))« 

Next^ in region 1 the viscous dissipation and heat-conduction terms 
are again neglected in the energy equation^ and terms that vanish by 
reasons of symmetry along the stagnation streamline are dropped. Thus 
the energy equation for region 1 takes the same form as equation (23) for 
the two-dimensional problem and; since the interface velocity is small; 
the interface temperature Tq is again approximately the stagnation tem- 
perature T-t . 

In region 2; the heat -conduction terns in the energy equation pre- 
dominate; and the equation reduces to the three-dimensional Laplace 
equation in the variable t) 

6x^ 6y^ 



In order to fit the boundary conditions on a spherical surface; the 
solution is given in teniis of spherical coordinates (f ; 0; and cp) • The 
general solution which preserves symmetry about the x axis (i.e.; which 
is independent of T) is 

00 

T = A + I + ^ (Cnr’^ + Pn(cos 0) (C2U) 

n=i 


where Pj^(cos 0) is the nth order Legendre polynomial in cos 0. If it 
is required that q be a constant, on the surface of the body, equa- 
tion (C24) can be reduced, on the stagnation streamline, to 


n = \ 




(C25) 
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then expanding in terms of e 


- 1 < < 1, results in 

Rb 


n = + L(e - e^) + . . . 


(C26) 


where L is 


00 

. B., y 

Rh L 


Cn(2n + 1)R^ 


n 


n=i 


Neglecting the q^uadratic term in 


€, evaluating t] at the interface, and transforming to the variable x, 
one obtains for equation (C26) on the stagnation streamline 


^ = 



X 

^b 


(C27) 


Then the rate of heat transfer to the stagnation point is 




Xb ^b 


(C28) 


which is identical in form with the zero-yaw solution for the two- 
dimensional-flow problem (eq. (31)) • Note, however, that the second 
derivative of pressure given by equation (B18) is larger by a factor 
of 4/3 than it is for the corresponding two-dimensional-flow case with 
the same shock-wave curvature (eq. (B17))*, Thus Xj^ given by equa- 
tion (C22) is changed by the factor (3/4) and the rate of heat trans- 
fer to an axially symmetric region becomes 





1/4 


1)J 


^ 5k 

.^^o Ri 


1/2 „ 1/2 

Re 

00 
R-K 



(C29) 


while the corresponding expression for Nusselt n^mlber is 


Nu 





1/2 


Re, 


1/2 


Rs/ ko(To - T^) 



(C30) 
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APPENDIX D 


EXAMIMTION OF ANALYTICAL RESULTS AND ASSUMPTIONS 


A nijmber of assiomptions have been made in the theoretical analysis, 
and it is desirable now to show that the solutions obtained are both 
realistic and consistent with these assumptions. In particular, it will 
be shown that the presumption of a constant second derivative of pressure 
normal to the stagnation streajnline yields solutions for the distance 
between shock wave and body which are reasonably close to observed values. 
Secondly, it will be demonstrated that the u velocity throughout 
region 2 is indeed small, as assigned in the analysis, if the Reynolds num- 
ber is large enough for continuum flow conditions. In addition, it will 
be shown that for region 2 the viscous -dissipation terms due to the 
u and V component velocity derivatives are small compared to the heat- 
conduction terms in the energy equation, again provided the Reynolds mmi- 
ber is not too small. These findings, then, help to justify the manner 
in which the momentum and energy equations were treated in the analysis. 

Now it is obvious that the assumption of an abrupt transition from 
nonviscous, convective flow to viscous, conductive flow is a substantial 
idealization of the actual flow.^ It is possible, however, to make a 
gross check on the self-consistency of this model by comparing the amount 
of heat convected across the interface with the amount conducted to the 
body surface. When this is done it is found that from a heat-flow point 
of view, the model is self-consistent (i.e., heat convected provides for 
heat conducted) . 

As a final point, a comparison will be made between the analysis of 
this paper and the heat-transfer solutions for low-velocity flow given by 
Howarth (ref. 7) and Cohen and Reshotko (ref. 8). 


Consider first axially symmetric flow. The velocity in region 1 was 
foimd to be (eq. (CIO)) 


^Strictly speaking, this idealized model should be considered simply 
a first approximation to the correct situation. A second approximation 
would be to divide the domain between the body and shock wave into three 
regions rather than two as was done here. 


Distance Between Shock Wave and Body 



(Dl) 
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Then, the shock-wave coordinate must he 


Xo = 


^s' - ^o' 


Ur 


(D2) 


It can be shown from the relations in Appendix B and equations (Cll) 
and (C12) that 


Us 


U, 


Ur 


^U„ 




(7g + 1)RS 


J P 3r^ 

8(3 - 27 ) U, 




(7s + l)Rs 


J2(7s - 1) > 


7.-1 R 


00 

2 


(D3) 


Substituting these relations into eqmtion (D2) yields 

X. ( 7 ^ - l)[l - J 2(7„ - 1)] 


R. 


2(3 - 27s) 


(d4) 


Note that for 7g = 1.5, Uq" vanishes and the velocity profile becomes 
linear. For this case Xg/Rg reduces to ( 7g - l)/4. 

The actual distance between the body and the shock wave is, of course, 
the sum of Xg and X'j). However, it can be shown from equations (C22) and 
(d 4) that x-^j is small compared to Xg for reasonably large Reynolds 
numbers, and xi^ will therefore be neglected. The ratio Xg/Rg calcu- 
lated from equation (d 4) for 7g equal 1.4 is O.IO 5 . Measurements of 
Xg/R-b taken from spark photographs of high-velocity spheres presented 
by Charters and Thomas (ref . 15) and Dugundji (ref . I 6 ) approach this 
value closely at high Mach nvimbers (i.e., Xg/R-^ about 0.11 at Mach nim- 
ber 4) . Heybey (ref. 17) has developed a theory which fits the data of 
references I 5 and 16 closely and, for the limit of infinite Mach number, 
predicts Xg/R-^ about 0.12. Thus it is seen that at high Mach n^mbers, 
the assumption that the second derivative of pressure is constant and 
that the ratio R^/Rs near \inity yields results which are consistent 
with experimentally observed distances between shock wave and body, as 
well as with the theory of Heybey. 
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It is of interest to calculate the shock-wave coordinate for two- 
dimensional flow as well. Recall that the solution for velocity in 
region 1 for this case is 


u 


Uo' 

sinh Cx 

C 


(D5) 


and thus 


— = ’ cosh Cx 

ax ° 


The velocity derivatives at the shock wave and at the interface, given 
in Appendix B, are, respectively. 


(7s + l)Rs 

> 

u ' = - / 1 ^ ^00 ^^(^S ~ 

° J P (7g + l)Rg ^ 

Then the product Cxg is given by 


(d6) 


CXs 


= arc cosh — === 

- 1 ) 


(D7) 


With CXg known and the velocity at the shock wave 


^s 


hJL 

7s + 1 


(d8) 


The shock-wave coordinate becomes 


Xc 




7s - 1 


sinh Cxc 


(D9) 
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For a 7g of lA, Cxg takes the value 0.75 since the sinh f\mc- 
tion is very nearly linear over this range ^ rather close hoimds on the 
shock-wave coordinate are in 5 )osed hy 


^S 

< 

^s 

< 

Us 



Us' 



^o' 


or 


^s 


2 


1 


< 




(DIO) 


(Dll) 


The exact theoretical solution for Xg/Rg at 7g = 1.4 is 0.236. 
According to the theory then, a shock wave with given radius of curvature 
should he detached from a cylindrical body about twice as far as from a 
sphere . 


Magnitude of Velocity in Region 2 


The y momentum equation in region 2 was reduced to 


M- 


6x 


6y^ 2 


+ l^o'^o 


(D12) 


The left side of this equation may be approximated by 1 (pu) with the 

OX 

presumption that velocity in region 2 is small. Then eq.uation (D12) may 
be integrated to 


p.u 


3 

o p X , 

6 


(D13) 


Solving for Uq, noting that velocity vanishes at x-j^, and substituting 
from equations h^) ( 19 ) } obtains 


u 


o 


2 

3 


^o'^b 


2 I 2PqUo' 

3 V ■ Ps 


(Dl4) 
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It follows that the ratio of interface velocity to the velocity at the 
shock wave is given by 



8 V ^ l6/t^oA Rb ^o' 

9 Ps^s ''s ~ 9 VPoo/Re„ 


(D15) 


which on substituting the relations given in Appendix B for u^' and Ug 
becomes 


/ 6 /do 1 

V^s/ 9 7 7s - 1 ^Rs/Reco 


(D16) 


It can be seen that for large Reynolds numbers^ of the order of 
hundreds or greater^ the velocity at the interface is small compared to 
the velocity at the shock -wave . Since the velocity in region 2 is every- 
where less than at the interface (see eq.. (D13))^ the solutions obtained 
for velocity are consistent with the assumption that velocity is small 
throughout region 2. 


Viscous Dissipation in Region 2 


Although the derivative of velocity vanishes at the body surface, it 
increases parabolically (see eq,. (D12)) to Uq’ at the interface. Since 
viscous dissipation terms due to this velocity shear were neglected in 
solving the energy equation, it will be shown that the maxim-um value of 
these terms, which occurs at the interface, is small compared to the heat- 
conduction terms like (note that by continuity Sv/Sy contributes 

a dissipation term of the same magnitude as ^u/^x) . From equation (30) 
it can be seen that the term S^q/Sx^ is nearly constant everywhere along 
the stagnation streatmline in region 2. Then the ratio of differential 
terns in the energy equation is, by equations (30) and (D14) , 


Vo(uq')^ 



(D17) 


If equation (Dl?) is evalijated for constant heat capacity and thermal 
conductivity proportional to the nth power of temperature , there is 
obtained 


Vq(^o')^ 


= l8(n + 1) 



(Dl8) 
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Substituting for velocity ratio Uq/u^o from equations (D16) and (b^) 
and for the ratio Rb/^b from equation (21) , there results 


(d^Ti/dx^)b 


= l6(n + l)Pr 






-1/2 

(DI9) 


Once again the square root of Reynolds niomber is the predominant term for 
conditions of continuum flov and thus the viscous dissipation terms in 

the energy equation are small compared to the conduction terms in region 2. 


Heat Convection Across the Interface^ x = 0 


Next consider the ratio of the heat convected across the interface. 
To 

CpdTj to the heat transfer at the stagnation point of the body, 

-q-jj'. The value of Uq given by equation (Dll4-) and -q, ' from equa- 
tion (32) yields 



P^o 



(D20) 


Again evaluating for constant heat capacity and the nth power tempera- 
ture function for thermal conductivity, and noting from equations (1^4-) 
and (19) that Uo'xb^ reduces to -(2no/p)j one obtains 


P^pCpTo 


= — Pr(n + 


1 ) 


(D21) 


This ratio is the order of unity, and thus the right magnitude of heat 
is convected across the interface to balance the heat conducted to the 
body. The above result also provides a check on the value of x-j^ which 
was obtained by matching vlq ' as a boundary condition of the y momen- 
tum equation. 
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Lov-Velocity Heat Transfer 


For hypersonic velocities it was found that taking shock-wave curva- 
ture equal to body curvature on the stagnation streamline gave approxi- 
mately the correct answer for the distance between the body and the shock 
wave^ so presumably the ratio should be taken near unity when cal- 

culating the heat transfer as well. Undoubtedly this ratio will be some- 
what less than unity for low Mach number supersonic flow^ and it is of 
interest to see what the solutions developed in this paper will predict 
for this case (even though the assumptions made in the analysis are not 
expected to hold as well for the low-velocity flow conditions). For this 
purpose it is convenient to express the body coordinate x-|^ in terms of 
(Sv/Sy)o which by continuity equals -Uq*. From equations ( 9 ), (l4), and 

( 19 ) 

then solving for Nusselt number from equations ( 32 ) and (3^) for the case 
of the cool wall (T^/Tq < < 1 , and n = l/2) one obtains 


2 Db , /p(Sv/^y) 

m = ^ — =0.U7]V / ° 

3 Xb Po 


(D23) 


The method of boundary-layer solution for low- velocity flow about a cylin- 
der given in reference yields for the derivative of velocity component 
normal to the stagnation streamline at the edge of the boundary layer, in 
the notation of this paper, 


^ 3-8U3 

\^yj i>b 

Substituting in equation (D23) results in 

^ 1/2 
Nu = 0.92 Reg 


{D2k) 


where the s tt mI 1 differences between |ig and |j.q are neglected. The con- 
stant 0.92 compares favorably with the value 0-95 given by Howarth for 
Pr = 0 . 72 . This agreement is especially remarkable in the light of the 
fact that the analysis of reference 7 is for constant thermal properties, 
while variation in thermal properties is an essential feature of this 
analysis . 
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Cohen and Reshotko (ref. 8) find that the solution for a compressible 
boundary layer gives the following relation at the stagnation point of an 
axially symmetric body 


^ = 0.440 

Db ^ 


(D25) 


for the case of a cool wall and a Prandtl nimber 0.?. If the radial com- 
ponent of velocity v is taken proportional to y, the ordinate can be 
eliminated and equation (D25) reduces to 


Nu = 0.440 Db 


p(8v/dy)Q 


(D26) 


The factor 0.44o given by Cohen and Reshotko compares favorably with the 
factor 0.47 given in equation (D23) . 
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APPENDIX E 

ESTIMATE OF RESERVOIR TEMPERATURES IN THE HYPERSONIC GUN TUNNEL 


A highly noniseutropic compression process, like the multiple reflec- 
tion of shock waves that occurs in the reservoir section of the hypersonic 
gun tunnel, is capable of producing high final air temperatures. An esti- 
mate of the magnitude of the reservoir temperatures attainable in the gun 
tunnel is obtained by equating the internal energy gained by the reservoir 
air to the work done in compressing this air. Thus assuming no heat loss 
(i.e., adiabatic compression) we have 



where the subscripts 1 and 2 refer to conditions before and after the com- 
pression, V is the volume of the gas, and Cy is the specific heat at 
constant volume. In terms of an average specific heat, Cy, and a mean 
pressure, p, equation (El) is 


- Tj = p(Vi-V2) 


(E2) 


By substitution of the mass conservation relation, ideal 

gas law p = pRT, and an average ratio of specific heats 7 defined by 
7 - 1 = R/Cy, equation (E2) can be transformed to 


^ _ 1 + (7 - l-)p/Pi 

1 + (7 - l)p/p 2 


Now the mean pressure during the compression process will probably be less 
than the final pressure, Pg, but greater than the average (p^ + ■p^)/2. If 
these values are taken as the limits of p, and in addition Pj^ is very 
small compared to P 2 ^ the temperature limits given by equation (E3) 
become approximately 


y - ^ < ^2 < ^ ^ £2 

7 + 1 Pi “ “ 7 Pi 


(e4) 
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For an average ratio of specific heats 7 = 1*3^ the limits become 


0.13 — < ^ < 0.23 — (E5) 

Pi “ Pi 

The compression ratio realized in the hypersonic gun tunnel is presently 
about 32 . Then for an initial temperature T^ = 530° R, the reservoir 
temperature limits indicated by equation (E5) are 


2200° R < T 2 < 3900 ° R 


(e6) 


Recall that the reservoir tempera tiure deduced from the velocity measure- 
ments is 2200° R^ the lower limit given by equation (e 6 ) . Undoubtedly 
the temperature corresponding to no heat loss from the reservoir would 
be somewhat higher than this, however. It is interesting to note that 
an isentropic compression of the same value, with no heat losses, would 
yield a final reservoir temperature Tg of only II 80 R. 
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APPENDIX F 

PROBLEMS IN CONNECTION WITH CALIBRATING 
THE HYPERSONIC GUN TUNNEL 

Static -Pres sure Measiirements 


The static pressvire system is considered to respond exponentially 
to the true static pressure vhich is in t\irn assiamed to be proportional 
to the pitot pressure. In this event we may write 


dp 

dt 



(FI) 


where t is the time constant of the system, and K is the ratio of 
static to pitot pressure. By integration of equation (F1) , the constants 
T and K stand in the linear relation 


/ P 
o 

p(t) - p(o) 


o 

p(t) - p(o) 


- T 


(F 2 ) 


The above integrals of pressure are evaluated as a function of time from 
the measured static-pressure and pitot -pres sure data. Then the values 
of K and t axe determined by the linear regression giving the least 
mean squares fit. It was found (see fig. I3) that a time constant 
T = 0.060 second and a ratio of static to pitot pressure K = 8.06x10“^ 
fit the pressiore data shown in figures 5(c) and 5(d) within a few percent 
over the entire interval of flow time from 0 to O.3 second. Accordingly, 
it is indicated that the value of p/H, corrected for time lag in the 
static pressTore measuring system, is 5 •32x10“® over this time interval. 


Condensed Phase Impurities 


It was observed during the course of experimenting with the hyper- 
sonic gun tunnel that the condensed phase ingsurities are absent from the 
shadow of models in the air stream. This observation suggests that the 
energy carried to a model by these inrpiirities may vary in proportion to 
the frontal area of the model. In this event, the heat transfer should 


CONFIDENTIAL 


mCA EM A55E02 


COWFIDENTIAL 


59 


tend to vary directly as the Reynolds number for the case of circular 
cylinders in the same stream. However, according to both theory and 
experiment (see fig. 9) heat transfer varies like the square root of 
Reynolds nimiber. There is added reason, then, to believe that condensed 
phase impurities have but a small influence on the initial heat-transfer 
rates presented in this paper. 


CONFIDENTIAL 


6o 


COKFIDENTIAL 


mCA RM A55E02 


APPENDIX G 

EFFECTS OF DISSOCIATION ON HEAT TMNSFER TO 
A CYLINDRICAL STAGNATION REGION 


Tempera t\ores in the disturbed flow about a vehicle in hypersonic 
flight may be sufficiently large (of the order of 5000° R or greater) to 
dissociate air molecules into atoms. Rates of dissociation are not pres- 
ently known with any degree of certainty, however, and it is therefore 
difficult to estimate whether the extent of dissociation will be negli- 
gible, partial, or in complete equilibrium in the time scale of flow 
about hypersonic vehicles of practical size . Since equilibrium disso- 
ciation is a limiting case, it is this situation which will be treated 
here . 


To a first approximation then, the temperature ratio across the 
shock wave is determined by the relation 


^s 


U 2 7 r(T )T M 2 

^00 _ ' OO ^ 00 00 00 

2 2 


(Gl) 


since the specific enthalpy, hg, is considered a known function of tem- 
perature for a given pressure. If the free-stream static pressiire is 
neglected, the momentum equation becomes 


P.. + = p U ^ 

-^s S S *^00 OO 

and substituting from the continuity relation 

p u = p U 

•^S S 00 00 

and the equation of state 

P = pR(T)T 


(G2) 


(G3) 


(g4) 
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the ratio of densities across the shock wave becomes 


/l.SiM (G5) 

Ps 2 7 2 


Then if 7 is defined by^ 

o 


p 

-22 = 


the pressiire ratio across the shock wave is 


- 1 


^s + 1 


(g6) 


Pc 


27 M 

00 oo 

7 +1 


(G7) 


From the equation of state (G^) and equation (g 6) , it follows that an 
equivalent expression is 


Ps _ ^ ^ 

p 7 - 1 T R(T ) 

where the ratio of specific gas constants R(Tg)/ R(T,jq) is also considered 
a known fimction of tempera tiure for a given pressure. Equations (GT) 
and (g 8) are solved for 7g ; then the rate of heat transfer to a cylindri- 
cal stagnation region may be calculated from equation (33) • Note that it 
is unnecessary to determine 7g with great precision, since it is the 
fourth root of (7g - l) which appears in equation (33) • 

The above calculations were made for a cylindrical stagnation region 
for the case of equilibrium dissociation at one atmosphere pressure and 


(g8) 


^This definition is used merely for convenience, so that the con- 
stants given in the heat-transfer solution (e'q. (33)) do not change in 
form. 
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free-stream temperature, Tio, equal 400° R.^ The enthalpy and specific 
gas constants were taken from data tabulated by Kreiger and White 
(ref . l8) and the transport properties were approximated by simple 
kinetic theory relations for a gas composed of hard spherical particles 
(see ref. 19). The quantitative accxiracy of this approximation is not 
high, of course, and the analysis cannot be expected to yield more than 
the order of magnitude for dissociation effects. The results of these 
calculations are presented in figure 12 both for flow in a state of 
equilibriimi dissociation and for flow in which dissociational energy is 
not excited. The calculated rates of heat transfer are higher when disso- 
ciation occurs but not in proportion to the large order of magnitude 
increase in transport properties. The explanation for this result is 
that the increase in transport properties is largely compensated for by 
a reduction in temperatiore level as the random kinetic energy is absorbed 
in the dissociational energy modes . In fact , as a resiilt of this com- 
pensation, Beckwith concludes on the basis of a boundary-layer analysis 
(ref . 20) , that the heat-transfer rate per unit area at the stagnation 
point is only slightly affected by equilibrium dissociation. A close 
comparison of the results of this paper with Beckwith's resialts woiild 
not seem warranted until the properties of dissociated air are more 
accurately known. 


^These conditions correspond roughly to upper atmosphere flight con- 
ditions (i.e., about 100,000 to 150,000 feet altitude at Mach n\jmbers from 
10 to 20) . At lower altitudes, higher pressures will be encoimtered 
behind the shock wave and the equilibrium dissociation will be less . It 
should be remembered, of course, that equilibrium dissociation may not 
occur in actual flight if the rate process is slow enough (see, e.g., 
ref. 21 for a discussion of dissociation rates) . 
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Nozzle valve opened — 



Figure /.- Schematic diagram of the hypersonic gun tunnel. 
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Figure 2- 20-mm hypersonic gun tunnel. 
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Valve seat 


(a) Photograph of nozzle valve and nozzle. 
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Sonic throat section 
.020 in. dia. 



Nozzle valve Nozzle 

(b) Schematic diagram. 


Figure 3.- Nozzle and nozzle valve for the 20-mm hypersonic gun tunnel. 
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Retaining yoke 



Strain gage elements 
housing 


Piston 


\ 

Insert to gun 
barrel 


A- 19913. 1 

(a) Reservoir pressure gage, gage insert, and yoke. 
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(b) Pitot pressure probe and capacitance pressure cell. 


Test station 

Pitot pressure probe 
Pressure cell 


Figure 4.— Calibration equipment for the 20-mm hypersonic gun tunneL 
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(c) Static pressure probes 



A- 19916. 1 


(d) Spark electrodes mounted in test section 


Figure 4- Concluded. 
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Figure 5- Reservoir, pitot and static pressure in the hypersonic 
gun tunnel as a function of time. 
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Figure 5— Concluded. 
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At = 44 ^ sec At - 5t^ sec 



(a) Schtieren photographs of the spark discharge 
disturbance sweeping through the test region. 



Figure 6.- Free-stream velocity as a function of time 

(At -time interval between spark discharge and Schiieren 
exposure, t-time of flow measured from nozzle valve 
opening). 
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040 in. dia. 


.020 in. dia. 


.012 In. dia. 



X)IO in. dia. 
.003 in. dia. 


A-19938.1 


(a) Thermocouple cylinders. 


A- 19917. 1 



(b) .020-inch-diameter thermocouple cylinder mounted 
in the test section at 45"* yaw. 


Figure 7.- Heat-transfer models. 
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Figure 6- Temperature of the .003- inch-diameter cylinder at zero 

yaw as a function of time. 
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Figure 9.- Comparison of average Nusseit numbers for cyiinders transverse 

to flow. 
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Figure !0- Effect of yaw on heat-transfer rates to the stagnation 

region of a cylinder. 
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Figure //.— Comparison of theory and experiment for the effect of 
yaw on the rate of heat transfer to a cylinder (M^=9.6, To =2 20CfR 

Tt/To ^ .24, Re^=t.26 x iO^ per foot). 
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Mach number, M^ 


Figure 12.- E ffect of equilibrium dissociation on 
the rate of heat transfer to the stagnation 
region of a cooi cylindrical body (p - / atm., 
To, - 400^R.) 
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Figure 13.- Variation of experimentally determined values of 

f,(t) with fzft). 
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